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It is a fundamental task in strongly correlated electron 
systems to examine interdependence among spin, charge, 
orbital, and lattice degrees of freedom. Recently, the ex- 
aminations of the interdependence under nonequilibrium 
conditions become quite urgent issues in various on-going 
subjects. Ultrafast transient optics is an approach ex- 
ploring new functionalities of materials and observing 
properties hidden in equilibrium conditions. The most 
important subject is photoinduced phase transition of 
low-dimensional transition metal oxides and organic ma- 
terials with strong electron correlation. 

One-dimensional (ID) Mott insulators, Sr2Cu03 and 
halogen-bridged Ni compounds, show photoinduced 
insulator-to-metal transition accompanied with its pi- 
cosecond recovery to the insulating state. ^'^ The time 
scale of the recovery is three order of magnitude faster 
than that for semiconductors.'^ The ID Mott insulators 
also exhibit gigantic third-order optical nonlinearity, and 
because of these two properties, they are promising fu- 
ture opto-electronics materials.^ The complete descrip- 
tion of optically excited Mott insulators is thus desired, 
but key theoretical concepts in nonequilibrium are still 
under construction. 

The photocarriers of the ID Mott insulators are called 
holon and doublon representing empty and doubly occu- 
pied sites, respectively. A holon and a doublon recombine 
by emitting energy to other elementary excitations. A 
problem is to clarify a pass way of energy dissipation due 
to the recombination. Two possible candidates are spin 
and phonon excitations, since antiferromagnetic (AF) ex- 
change energy and highest phonon frequencies are of the 
same order (■~0.1 eV).^ Since high-energy states created 
by optical excitation may violate the separation of spin 
and charge degrees of freedom inherent in ID electron 
systems,^ a pass way for energy dissipation through a 
spin channel can be expected. However, recent numer- 
ical studies have shown robustness of the spin-charge 
separation for nonequilibrium steady states.^'* It is thus 
necessary to make clear a coupling of spin and charge de- 
grees of freedom under photoirradiation. As for phonon 
relaxation, pump-probe experiments have been done for 
various TTF-TCNQ salts with different magnitudes of 
electron-phonon (EP) coupling.^ K- and Rb-TCNQ show 



spin-Peierls (SP) transition at T^ =395 K and 220 K, re- 
spectively, and their photocarriers are once localized as 
polarons at around 70 fs, and then recombine with a 
few ps. On the other hand, ET-F2TCNQ does not show 
SP transition, and metallic photocarriers decay within 
200 fs. Therefore, a fundamental question to be answered 
is about what is driving force of ultrafast relaxation of 
the ID Mott insulators when the charge carriers couple 
weakly with spin and lattice. Since EP coupling is also 
present in semiconductors, we need to answer another 
question why phonon relaxation in the Mott insulators 
is much faster than that in the semiconductors. 

In this Letter, we incorporate time dependent vec- 
tor potential of laser pulse into density-matrix renor- 
malization group (DMRG) simulation of a ID Hubbard- 
Holstein model to answer the questions raised above. Ex- 
amining initial relaxation after irradiation, we find that 
the spin-charge coupling exists for strong excitation in 
the case without phonon degrees of freedom. With in- 
troducing even small magnitude of EP coupling, the ef- 
fect of the spin-charge coupling is suppressed, and many 
phonons are excited dynamically in the system. Phonon 
bases required for time-dependent calculation are much 
larger than those for the calculation of the optical con- 
ductivity with the same magnitude of EP coupling. This 
dynamical generation of phonons characterizes phonon 
relaxation in Mott insulators. We also discuss the differ- 
ence of relaxation dynamics between Mott insulators and 
semiconductors . 

We start with a ID extended Hubbard-Holstein model 
with the classical vector potential of pump light: 
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where c] ^ (q,ct) is a creation (annihilation) operator of 
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an electron at site i with spin a, and bj {hi) is a creation 
(annihilation) operator of a phonon at site i. This model 
includes electron hopping, t, on-site and nearest-neighbor 
Coulomb repulsions, U and V , respectively, phonon fre- 
quency, ujQ, and EP coupling, g. The time-dependent vec- 
tor potential, A(t), is defined by 



A{t) = Aoe-(^-^")'/(2^-') cos(c.pu„,p(T - To)), 



(2) 



where t^ is the duration time and Wpump is the frequency 
of pump laser. We note that, in the present study, both of 
electronic and phononic degrees of freedom are quantized 
in contrast to a recent work where phonons are treated 
as classical lattice vibrations. ^'^ 

We calculate linear optical absorption spectrum before 
pumping, which is related to the current-current corre- 
lation function given by 
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with the current operator j = itY^^ ^{c[ ^Cj^^ „ ~ H.c. 
and the ground state |0) with energy Eq. 

Photodoping rate due to pumping, (5, which is con- 
trolled by the amplitude of laser pulse Aq in (2), is de- 
fined by 



5 = lim {E{t) - -Eo)/(La;pu,np) 



(4) 



where L is the system size and E{t) — {t\H\t). \t) is 
the solution of the time-dependent Schrodinger equation, 
i{d/dT) \t) ~ H{t) |r), and the solution is given by 
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with T being the time ordering operator. 

We calculate the time evolution of doublon number, 
N(i{t), phonon number, 7Vph(T), and local spin correla- 
tion, Cs(t). They are defined by 
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For numerical calculation of these time-dependent 
quantities, we apply the DMRG technique as well as the 
numerically exact diagonalization method. In the DMRG 
technique, we perform a DMRG run for a given time r 
and repeat this for many r. Here, we optimize just one 
time-dependent wave function |r) for a given time, and 
thus can reduce the number of the target states in the 
reduced density matrix, leading to better performance. 
For the calcination of |r), we recursively use the following 
equations for small time slice Ar, 
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Fig. 1. Current-current correlation function x('^) of ^ H^ ex- 
tended Hubbard-Holstein chain for L = 12 and a broadening 
factor 7 = 0.2i. The fine solid, bold solid, and bold dashed lines 
correspond to g = 0, O.li, and 0.2i, respectively. 



until we obtain |r). Here, M represents the number of 
the Lanczos step, t„ = (n — 1/2)At, \j) and Ej are the 
eigenvector and eigenenergy of the tri-diagonal matrix 
for H{Tn), respectively, and {j\ (n — 1)At) is the first 
component of \j). Note that it is more sophisticated to 
do continuous time sweeping based on the adaptive time- 
dependent DMRG method. ^^ However, the wave function 
transformation crucial for this method only leads to poor 
numerical precision in the present model, and we do not 
use the adaptive method. 

We consider half-filled systems, and the maximum sys- 
tem size is L = 12. In our DMRG calculations, electron 
and phonon degrees of freedom are represented as differ- 
ent sites and an open boundary condition is used. The 
parameters U, V, and wq in (1) are taken to be lOi, 
2t, and O.lt, respectively.^^ Among parameters in (2), Td 
and To are fixed to be t^ — 2/t and tq = 5/t, respectively. 
The truncation number m of density-matrix eigenvalues 
in the DMRG procedure is taken to be up to to = 1000. 
The number of phonon is taken up to 6 states per phonon 
site. 

We first show x(cj) for various g values in Fig. 1. The 
calculation was done by the dynamical DMRG and the 
obtained results agree qualitatively with previous numer- 
ical results. ^^ Because of small system size, xi'^) exhibits 
discrete peaks with finite energy difference. For g = O.li, 
line shape with a broadening factor 7 = 0.2Hs similar to 
that without g. In this case, it is enough to prepare two 
phonons on each site. In the following, we take g ~ O.lt 
and g = 0.05t in order to make clear whether weak EP 
coupling affects the time evolution very strongly. Here, 
the effective EP coupling is defined by A = g'^/{2tujo), 
and the value of A is A = 0.05 and A = 0.0125 for g = O.lt 
and g — 0.05i, respectively. 

Before going to the extended Hubbard-Holstein model, 
we study the time evolution of the extended Hubbard 
model (g — 0) to make clear how the coupling between 
photocarriers and spins is induced by strong laser pulse. 
In Fig. 2, we show time evolution of E{t), N^ir), and 
Cs(t). The photodoping rate S is set to be 0.37 corre- 
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Fig. 2. (a-c) Time evolution of (a) E{t), (b) N^{t), and (c) Cs(t) 
in a L = 10 extended Hubbard chain. Pump-laser irradiation 
continues up to r ~ 10/t and photodoping rate 5 = 0.37. In- 
sets show enlarged scale of the corresponding quantities. Here, 
roughly 4 photons come in and then 2 doublons are created. Be- 
cause of high density excitation, the number of photon and the 
number of doublon-holon pair are not the same, (d) System-size 
dependence of A^d(''') by exact diagonalization. Bold solid line 
for L = 10 and fine solid line for L = 8. We take Wpump = 7t 
(close to the energy of the excitonic bound state) and Aq = 1.0 
(<5 ~ 0.26). 



spending to strong excitation. This case roughly corre- 
sponds to an experimental situation where clear Drude 
weight appears in absorption spectra just after pump- 
ing. Wpump is set to he U — V = 8i, and then the opti- 
cal pump mainly excite free carrier states. E{t), iVd(T) 
and Cs(t) gradually increase during photoirradiation as 
expected. After the pumped pulse has been turned off 
(r > 10 /t = 2to), E{t) is hardly changed because of 
treating a closed system. 

We find in Figs. 2(b) and 2(c) that N(i{T) decreases 
as Cs(t) increases with time (see insets) accompanied by 
oscillating behaviors. Unfortunately, the time resolution 
is limited in our DMRG due to the sparseness of time 
mesh. Thus, in order to see clearly the meaning of the 
decrease of N^ir) and the oscillation, we show the exact 
diagonalization results of A^dC^) for L = 8 and L = 10 in 
Fig. 2(d). We take a pameter set different from that for 
Fig. 2(b), so that the oscillation looks clearer. Note that 
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Fig. 3. Time evolution of various quantities in an extended 
Hubbard-Holstein chain with L = 12. The bold solid, the fine 
solid, and the dashed lines represent the cases with g = O.lt, 
g = 0.05t, and 5 = 0, respectively, (a) E{t) - Eq, (b) Afph(-r), 
(c) N^(t), and (d) Cs(t). The inset of (b) represents the first 
derivative of A'^ph (''')■ We take Wpump = 6.75t corresponding to 
the excitonic bound state (see Fig. 1), and Aq = 2.0. Pump-laser 
irradiation continues up to r ^ 10/t. 



photoexcitation creates bound excitons in this case, and 
due to this fact the recombination rate is larger than that 
for Fig. 2(b). We see that the osciUation depends on the 
system size, while the decrease of N^ir) does not. There- 
fore, we can safely say that the decrease of N^ir) cap- 
tures an essential feature of relaxation. The oscillation is 
composed of two periods. One is Wp^jj which is much 
faster than our time mesh. The other one is roughly L/t 
which can be observed in the DMRG data. Physically 
this time scale corresponds to electron motion from one 
edge of the chain to another edge. This is also related to 
descrete energy level structure of finite-size system. 

We have confirmed that the gradual increase of Cs(r) 
seen in Fig. 2(c) is not clearly obse rved in cases of weak 
excitation.^ The decrease of N^ir) and the increase of 
C's{t) mean that holon and doublon created by pump 
pulse recombine with each other and the recombination is 
accompanied by an energy flow from the photocarricrs to 
the spin degree of freedom. However, it should be noted 
that the slope of time dependence in Cs(t) is very grad- 
ual. Furthermore, the change of Cs(r) from tr — 10 to 50 
(^0.2) is smaller than the value of one spin flip (~1.0). 
Therefore, relaxation through the spin channel is not ef- 
ficient in the extended Hubbard model, even though the 
spin-charge coupling becomes evident as the pump power 
increases. This is qualitatively consistent with recent re- 
ports showing robustness of the spin-charge separation 
for nonequilibrium steady states.^'® 

The effect of EP coupling on relaxation dynamics is 
examined in Fig. 3, where CLipump is tuned to the ex- 
citonic bound states shown in Fig. 1. Until the pump 
pulse is tuned off, all of charge, spin, phonon degrees 
of freedom are disturbed. After pumping (tr > 10), the 
doublon number N(i{t) decreases as is the case without 
EP coupling. In the latter case, the resulting energy loss 
in the charge degree of freedom was transferred to the 
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spin channel accompanied by the decrease of spin corre- 
lation, i.e., increase of effective spin energy. In the former 
case with EP coupling, spin correlation is also reduced. 
However, the reduction of spin correlation is suppressed 
by increasing EP coupling, as is seen in Fig. 3(d). This 
implies the suppression of spin-charge coupling due to 
the presence of EP coupling. 

The number of phonons increases after pumping as 
shown in Fig. 3(b). A^ph(T) does not contain oscillation, 
since the inverse phonon frequency uJq = 10/t is slower 
than the other time scales. We find for g = O.lt that 
^ph(''') in Fig. 3(b) has two regions in the time domain 
after pumping {tr > 10): A'ph(T) increases rapidly and 
the slope changes at tr ~ 30 as is seen in the first deriva- 
tive of 7Vph(r) [the inset of Fig. 3(b)], followed by a gen- 
tle slope for tr > 30. This characteristic time of kink 
T ^ 30/t roughly corresponds to the duration time of 
laser pump (~ 10/t) plus the inverse of effective EP 
coupling, (At)^^ ~ 20/i. In a case with g — 0.05i, the 
inverse EP coupling is (At)^^ — 80/t, and the change 
in the slope is expected to occur at around r ~ 90/i. 
However, the kink structure is not clearly seen, since the 
coupling is too small. The rapid increase of phonon dur- 
ing 10/t < T < 30/t indicates an energy transfer from 
charge degree of freedom to phonon one in the initial 
stage of relaxation. In other words, initial relaxation is 
dominated by phonons. 

In the following, let us examine the mechanism of 
phonon relaxation in light of the presence of the two time 
regions. As shown in Fig. 3(b), the expectation value of 
phonon number per site is less than one, but the size 
of local phonon Hilbert space taken in the calculations 
is huge. Then, we have two possibilities for the relax- 
ation processes. The first one is that the recombination 
of a holon and a doublon is acceralated by the phonon 
emission, if these carriers were still mobile against the 
EP coupling. In this case, the relaxation has been al- 
most finished at around tr = 30. On the other hand, 
if the phonon generation makes these carriers polaronic, 
the time scale tr = 30 represents the polaron formation, 
eventually leading to slow recombination dynamics in the 
later time region. In the present case with g — O.lt, the 
former scenario is realized. We have tried preliminary cal- 
culation of a correlation function (t| bjbjni^-fni^i \t) in a 
6-site chain, and have actually confirmed that the spatial 
distribution of phonons around photocarriers is uniform. 
When we take g = 0.2t, the distribution starts to concen- 
trate on the photocarrier site. The detailed results will 
be shown elsewhere. 

Let us also estimate the typical value of the time 
ir = 30. Our parameter set roughly corresponds to that 
for ET-F2TCNQ, since g hardly affects the optical con- 
ductivity. The band width of the optical conductivity of 
our single-band model is 8t, and experimentally the main 
band is located at around 0.5 ^ 2.1 eV. Then, the elec- 
tron hopping is estimated to be 0.2 eV. In this case, we 
obtain r = 30/t ^ 100 fs. This value is not contradictory 
to the value 200 fs for ET-F2TCNQ. 

In our previous study for the effect of tZ-depcndence 
on charge-phonon coupling, we have found that photo- 



carriers are dressed with phonon cloud strongly as U in- 
creases for the large- 1/ region. ^^ The result is quite sug- 
gestive for the difference of phonon relaxation dynam- 
ics between band and Mott insulators. In the largc-C/ 
limit, the model can be mapped onto a holon-doublon 
model. In this model, one holon-doublon pair exists and 
the spin degree of freedom is completely traced out. In 
the sense that only the charge degree of freedom remains, 
a phonon effect on these particles can be matched to 
that on charge carriers photo-doped into a band insu- 
lator. Taking this correspondence into account, we can 
judge that the photocarriers in the band insulator couple 
with phonons more strongly than in the Mott insulators. 
Strong EP coupling would accelerate localized polaron 
formation, leading to slow relaxation dynamics. The dif- 
ference between band and Mott insulators is also seen in 
the DMRG calculation of spectral function, ^^ where we 
have needed a large number of local phononic states for 
band insulators in comparison with Mott insulators with 
the same EP coupling constant. 

Summarizing, we have studied the relaxation dynamics 
of photocarriers in the ID Mott insulators with EP cou- 
pling. The EP coupling dominates the spin-charge cou- 
pling in the initial relaxation, even if the optical conduc- 
tivity is not affected by the EP coupling. We discussed 
the difference of relaxation between Mott and band insu- 
lators combining the present results with t/-dependence 
on polaron formation. 
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